Cavitating Langmuir Turbulence in the Terrestrial Aurora 
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Langmuir cavitons have been artificially produced in the earth's ionosphere, but evidence of 
naturally-occurring cavitation has been elusive. By measuring and modeling the spectra of electro- 
static plasma modes, we show that natural cavitating, or strong, Langmuir turbulence does occur in 
the ionosphere, via a process in which a beam of auroral electrons drives Langmuir waves, which in 
turn produce cascading Langmuir and ion- acoustic excitations and cavitating Langmuir turbulence. 
The data presented here are the first direct evidence of cavitating Langmuir turbulence occurring 
naturally in any space or astrophysical plasma. 

PACS numbers: 94.05. Lk, 94.05.Pt, 94.05.Fg, 94.20.wj 



Langmuir turbulence is known to occur in controlled 
laboratory [1, 2 and space plasma experiments [3-5 and 
is thought to occur naturally in a variety of space and 
astrophysical plasmas, including pulsar magnetospheres 
|6^, the solar corona [7^, the interplanetary medium [8 , 
planetary foreshocks [9], the terrestrial magnetosphere 
[ro, and the ionosphere plT[T3] . In its most developed 
form, this turbulence contains electron Langmuir modes 
trapped in dynamic density depressions known as cavi- 
tons ^14Hl6j . Cavitons have been shown to be artificially 
produced in the earth's ionosphere during high-power 
radiowave pumping experiments as deduced from radar 
spectra containing simultaneously-excited up and down- 
shifted Langmuir and ion-acoustic lines plus a central 
peak due to cavitation [3-5 , but evidence of naturally 
occurring cavitation has until now been elusive. 

Between 18 and 21 UT on 11 and 12 November 1999 
a measurement program designed to detect both ion- 
acoustic and Langmuir modes was run on the European 
Incoherent Scatter Scientific Association (EISCAT) 224- 
MHz radar located near Troms0 in northern Norway 
(local standard time in Norway is UT plus one hour). 
The principal objectives were to observe enhanced waves 
stimulated by high-power radiowave pumping, and, in 
the event of auroral activity, to gather data on nat- 
ural energetic waves [17 . On both nights conditions 
were disturbed, and enhanced echoes were detected, the 
strongest being on 11 November between 18:18:30 and 
18:21:30 UT, during the passage of an aurora through 
the vertically-directed radar beam. Fig. [l] presents pa- 
rameters derived from the ion-acoustic backscatter be- 
tween 18:15 and 18:28 UT, during the most intense au- 
roral event. Fig. |2] shows the intensities of Langmuir and 
ion- acoustic backscatter as a function of height and time. 
The prominent features occurring between 18:18:30 and 
18:20:30 UT and at 18:23:30 UT near 300 and 250 km al- 
titude, respectively, are backscatter associated with the 



aurora, and are the most energetic natural events ob- 
served on either night. Two other events occurred later 
that evening and two more on the following evening. 
Weak ion-acoustic enhancements occurred during each 
event; the Langmuir enhancements, however, are always 
stronger. Fig.|3]shows up and down-shifted spectral lines, 
or "shoulders", which are produced by Doppler-shifted 
backscatter from the down and up-going ion-acoustic 
waves, respectively. The shoulders are strongly enhanced 
indicating that the waves are nonlinearly amplified. In 
addition, there is a strong central peak, a feature not 
present in thermal-level spectra. 

The results of a computation made for plasma param- 
eters matching those which occurred during this obser- 
vation are shown in Fig. |4j A numerical code incorporat- 
ing a one-dimensional periodic version of the Zakharov 
equations was used |20], capable of producing the 
full range of cascading (sometimes called weak), coex- 
istence (transitional), and cavitating (strong) Langmuir 
turbulence. Energy was supplied by a downward- going 
flux, or beam, of electrons [21], which excites a Langmuir 
"pump" wave via the bump-on-tail instability. In the 
cascading turbulence scenario [22], the Langmuir wave 
then undergoes parametric decay into daughter Langmuir 
and ion- acoustic waves. These waves however exist only 
within two relatively narrow bands of wave numbers: the 
Langmuir band defined by the driving beam (see caption 
to Fig. 4) and the ion-acoustic band at about twice that 
value. Furthermore, a radar sees only the wave number 
that matches the Bragg scattering condition for the radar 
wavelength. This means that, for beam-driven cascading 
turbulence, a radar will see either enhanced Langmuir 
waves or enhanced ion- acoustic waves, but both may be 
seen simultaneously only when the velocity spread of the 
beam is sufficiently broad, approaching the absolute ve- 
locity of the beam itself. Cavitating turbulence is dif- 
ferent in that enhanced wave modes cover a range of k 
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FIG. 1. (color online) Background ionospheric parameters 
measured during the most prominent auroral event. The 
panels show, from top to bottom, electron density, electron 
temperature, ion temperature, and vertical ion velocity (pos- 
itive indicates motion away from the observer) at 1-min time 
resolution. Typical auroral plasma signatures can be seen, 
namely a sharp increase in electron density localized in time 
and space, a corresponding increase in electron temperature, 
localized small increases in ion temperature, and, in the ve- 
locity plot, high- altitude ion outflow to space [T9] . 



space which extends broadly on both sides of the wave 
number of the pump wave irrespective of the beam veloc- 
ity breadth. In the coexistence or transitional case the 
wave number spectrum extends below the pump wave 
number to zero, but dies out rapidly for Langmuir wave 
numbers greater than the Langmuir pump and for ion- 
acoustic wave numbers greater than twice the Langmuir 
pump. 

In the simulations presented in Fig. [4] several different 
beam energies are modeled, each capable of producing 
cavitating turbulence, and the spectra are those which 
would be seen by a radar of the same wavelength as was 
used for the observations; however, the beam parameters 
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FIG. 2. Incoherent scatter intensity profiles from up and 
down-going Langmuir and ion-acoustic waves recorded during 
the auroral event discussed in Fig. 1. Four distinct sources 
of scattering can be identified. (1) The dark background 
in the ion-acoustic channel is backscatter from thermal-level 
waves. (2) The relatively faint bands in the down-going and, 
more prominently, up-going Langmuir channels correspond 
to 3 and 5 MHz Langmuir waves weakly enhanced by low en- 
ergy, direct and backscattered diffuse electron precipitation 
associated with the aurora. (3) The repeated 10-s-long en- 
hancements seen in all three channels at about 225 and 125 
km before and after 18:18:30, respectively, are backscatter 
from waves enhanced by experimental 4.04-MHz high-power 
radiowave transmissions [T7]. (4) The intense features occur- 
ring between 18:18:30 and 18:20:30 UT near 300 km and at 
18:23:30 UT near 250 km are backscatter associated with the 
aurora. The top edges of the Langmuir enhancements give the 
approximate altitudes of the enhanced backscatter. Uncoded 
420-/is and 25-/is pulses were used to measure the Langmuir 
and ion-acoustic backscatter, respectively. 



were chosen so that the Langmuir and ion-acoustic wave 
numbers for cascading turbulence would not match that 
of the radar. The coincident enhancement in space and 
time of both ion-acoustic and Langmuir backscatter at a 
single radar wave number is a prediction characteristic of 
cavitating Langmuir turbulence, and constitutes critical 
evidence for its occurrence [13l |T4] . 

A second key feature of this observation lies in the 
shape of the ion-acoustic spectra, which consists of en- 
hanced up and downshifted "shoulders" and an enhanced 
central peak, shown in Fig. [3j The observed spectra 
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FIG. 3. Power spectral densities of the naturally-enhanced 
ion-acoustic backscatter showing enhanced shoulders and en- 
hanced central peaks at 2-kHz frequency and 3.75-km range 
resolution. A 475-/is phase-coded pulse with a 25-/is baud 
was used. The change in enhancement height with time is a 
result of the vertical pointing of the radar, while the auroral 
structure is oriented along the geomagnetic field, tilted 13° 
south of vertical, and drifting south; this effect can also be 
seen in Fig. 2. The progression from higher to lower heights 
with time corresponds to a drift velocity of about 80 m/s. 



match the computed spectra in Fig. |4] very well. En- 
hanced shoulders are a standard feature of all past ob- 
servations of naturally-enhanced ion-acoustic backscat- 
ter [23], but only one past result shows an enhanced 
ion-acoustic spectrum with a central peak [24 . The ap- 
pearance of the non-Doppler-shifted central peak indi- 
cates the presence of meter-scale, non-propagating den- 
sity wells known as cavitons. A central peak is not a 
necessary feature of cavitating turbulence but will oc- 
cur when the spacing between the cavitons matches the 
Bragg condition of the radar. The spacing, in turn, is 
roughly proportional to the inverse square root of the 
energy density of the primary Langmuir waves, or, in 
turn, of the pump beam [131 123 • This is not likely to 
be seen in all observations, both because of the matching 
requirement and because a relatively high beam energy 
density is required, which is most likely to be observ- 
able using long wavelengths such as in the observations 
reported here and previously |24j. 

Many other features of the observations may also be 
explained by the cavitating turbulence model. (1) Six 
events were observed, four on the 11th and two on the 
12th, and both Langmuir and ion- acoustic enhancements 
exist in all cases, but in all events the ion-acoustic en- 
hancement is weaker. This agrees with the computed 
examples. (2) The measured backscatter intensity in 
the down-going Langmuir channels is weaker than the 
up-going. Similar differences can be seen in the com- 
puted spectra. For the case of 18 eV, the radar sees 
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FIG. 4. Results from a 1-d simulation made for plasma pa- 
rameters matching those during this observation and with a 
downward-going beam of electrons at beam energies of 18, 42, 
75, and 170 eV. The beam creates a bump-on-tail distribution 
which excites Langmuir waves according to the resonance con- 
dition Vh ~ Al/l, where Vh is the beam velocity and Al and 
/l are the Langmuir wavelength and frequency. The spectra 
of up and down-going Langmuir and ion-acoustic waves were 
calculated for a wavelength of 0.67 m, matching the Bragg 
backscatter condition of the radar. Both the precipitating 
electron beam and the thermal background were included in 
the driving terms. The beam velocity spread ratio Avh/vh is 
0.3 in all cases, where Avh is the velocity spread. The beam 
density ratio nb/rip is 2 • 10~^, where rih is the beam density 
and rip is the density of the surrounding plasma, 3 • 10^^ m~^, 
as determined from the 5-MHz frequency of the enhanced 
Langmuir modes. The electron collision frequency is 200 s~^. 
These parameters are well within the ranges known to occur 
in the auroral ionosphere 21 . The peak value of each spec- 
trum is given at the top right of each panel. The Langmuir 
spectra are shown on identical log scales with arbitrary but 
equal reference levels and with a minimum value 10~^ that of 
the maximum; the ion-acoustic spectra are shown on a linear 
scale in arbitrary units with a minimum value of zero. 



down-going Langmuir modes at a wave number corre- 
sponding to the negative, or damping, slope of the ve- 
locity distribution function of the down-going beam, and 
heavy damping is seen. For wave numbers on the pos- 
itive slope of the beam the opposite will occur. (3) In 
five of our cases the ion-acoustic enhancement disappears 
before, or at the same time as, the Langmuir enhance- 
ments (within the 10-s time resolution of the observa- 
tions). However, in our strongest event, at 18:18:30 on 
the 11th, enhanced ion-acoustic backscatter is seen be- 
fore and after the enhanced Langmuir backscatter. In 
this case it is possible that the turbulence develops from 
coexistence to cavitating and back again as the driv- 
ing beam grows and then decays. In the coexistence 
regime Langmuir mode wavenumbers are cut off at a 
value roughly half that of the ion-acoustic mode, so the 
radar may see only the ion-acoustic enhancement until 
the cavitating turbulence is fully developed. A strong 
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event would also cause greater electron heating, reduc- 
ing damping of ion-acoustic modes and contributing to 
a longer ion-acoustic enhancement. (4) The relationship 
between beam energy and ion-acoustic damping can also 
account for the relatively strong ion-acoustic enhance- 
ment after 18:20:00 UT: the electrons had been heated 
during the event and the ion-acoustic damping rate re- 
duced, allowing the ion-acoustic backscatter to remain 
strong even as the drive began to weaken. 

A significant feature in the computed Langmuir spec- 
trum at 170 eV, and very weakly at 42 and 75 eV, is 
the appearance of a broad spectrum at and near the cold 
plasma frequency. This feature, which is not resolvable in 
the radar measurements presented here, is due to Lang- 
muir waves trapped in cavitons. It appears after a suffi- 
cient period of pumping at a sufficiently high level. The 
central peak and the broad spectrum, both due to cav- 
itation, may be seen under somewhat different circum- 
stances: The central peak requires a beam energy den- 
sity that allows the caviton spacing to match the radar 
Bragg condition, while the broad spectrum requires a 
fixed electron density with a beam velocity high enough 
and beam duration long enough to allow formation of 
trapped Langmuir waves matching the radar Bragg con- 
dition. Both are seen in high-power radiowave experi- 
ments [5 , in which a fixed pump frequency substitutes 
for a fixed electron density. 

The data presented here provide the first direct evi- 
dence of naturally occurring cavitating Langmuir turbu- 
lence, thought to be important in space and astrophysi- 
cal plasmas as varied as pulsar magnetospheres and the 
Earth's ionosphere [6-14 . Further observations of Lang- 
muir turbulence in the ionosphere may yield advances 
in our understanding of supra-thermal electron distribu- 
tions [26], naturally-enhanced ion-acoustic waves [El [22], 
natural ionospheric radio emissions [27] , anomalous resis- 
tivity [28l [29] , and auroral currents and dark aurora [18] . 
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